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We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross- 
section from future ton-scale direct detection experiments using germanium, xenon or argon as 
targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter 
experiments which will come online within 5 to 10 years, the degree of complementarity between 
different targets is quantified. We investigate how the uncertainty in the astrophysical parameters 
controlling the local Dark Matter density and velocity distribution affects the reconstruction. For a 
50 GeV WIMP, astrophysical uncertainties degrade the accuracy in the mass reconstruction by up 
to a factor of ~ 4 for xenon and germanium, compared to the case when astrophysical quantities are 
fixed. Ifowever, combination of argon, germanium and xenon data increases the constraining power 
by a factor of ~ 2 compared to germanium or xenon alone. We show that future direct detection 
experiments can achieve self-calibration of some astrophysical parameters, and they will be able to 
constrain the WIMP mass with only very weak external astrophysical constraints. 



I. INTRODUCTION 

Many experiments are currently searching for Dark 
Matter (DM) in the form of Weakly Interacting Mas- 
sive Particles (WIMPs), by looking for rare scattering 
events off nuclei in the detectors, and many others are 
planned for the next decade [IHB]. This direct DM de- 
tection strategy has brought over the last year several 
interesting observations and upper limits. The results of 
the DAMA/LIBRA [Tj and, more recently, the CoGeNT 
[8] collaborations have been tentatively interpreted as 
due to DM particles. It appears however that these re- 
sults cannot be fully reconciled with other experimen- 
tal findings, in particular with the null searches from 
XENONIOO ISHlI] or GDMS [T^, and are also in tension 
with ZEPLIN-III [T3j. In this context, the next genera- 
tion of low-background, underground detectors is eagerly 
awaited and will hopefully confirm or rule out a DM in- 
terpretation. 

If convincing evidence is obtained for DM particles 
with direct detection experiments, the obvious next step 
will be to attempt a reconstruction of the physical param- 
eters of the DM particle, namely its mass and scattering 
cross-section (see e.g. Refs. [HHIS])- This is a non-trivial 
task, hindered by the different uncertainties associated 
with the computation of WIMP-induced recoil spectra. 
In particular, Galactic model uncertainties - i.e. uncer- 
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tainties pertaining to the density and velocity distribu- 
tion of WIMPs in our neighbourhood - play a crucial role. 
In attempting reconstruction, the simplest assumption to 
make is a fixed local DM density po ~ 0.3 GeV/cm'^ and a 
"standard halo model" , i.e. an isotropic isothermal sphere 
density profile and a Maxwell-Boltzmann distribution of 
velocities with a given galactic escape velocity Vesc and 
one-dimensional dispersion cr^ = Vq/2 = {vq be- 

ing the most probable velocity and visr the local circular 
velocity, see below). However, the Galactic model param- 
eters are only estimated to varying degrees of accuracy, 
so that the true local population of DM likely deviates 
from the highly idealised standard halo model. 

Several attempts have been made to improve on the 
standard approach [TBHTO] . In the case of a detected 
signal at one experiment, recent analyses have studied 
how complementary detectors can extract dark matter 
properties, independent of our knowledge of the Galac- 
tic model PD]. Certain properties of dark matter may 
also be extracted under assumptions about the nature 
of the nuclear recoil events [5T]. Furthermore, eventual 
multiple signals at different targets have been shown to 
be useful in constraining both dark matter and astro- 
physical properties [Hj and in extracting spin-dependent 
and spin-independent couplings |231 124 ) . Here, using a 
Bayesian approach, we study how uncertainties on Galac- 
tic model parameters affect the determination of the DM 
mass and spin-independent WIMP-proton scattering 
cross-section cr^j. In particular we focus on realistic ex- 
perimental capabilities for the future generation of ton- 
scale detectors - to be reached within the next 10 years 
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- with noble liquids (argon, xenon) and cryogenic (ger- 
manium) technologies. 

The main focus of this paper is the complementarity 
between different detection targets. It is well-known (see 
e.g. [2]) that different targets are sensitive to different 



directions in the 



'SI 



plane, which is very useful 



to achieve improved reconstruction capabilities - or more 
stringent bounds in the case of null results. This problem 
has often been addressed without taking proper account 
of Galactic model uncertainties. Using xenon (Xe) , argon 
(Ar) and germanium (Ge) as case-studies, we ascertain 
to what extent unknowns in Galactic model parameters 
limit target complementarity. A thorough understand- 
ing of complementarity will be crucial in the near future 
since it provides us with a sound handle to compare ex- 
periments and, if needed, decide upon the best target to 
bet on future detectors. Our results also have important 
consequences for the combination of collider observables 
and direct detection results (for a recent work see |25|). 

Besides degrading the extraction of physical proper- 
ties like and cr^j, uncertainties in the Galactic model 
will challenge our ability to distinguish between different 
particle physics frameworks in case of a positive signal. 
Other relevant unknowns are hadronic uncertainties, re- 
lated essentially to the content of nucleons [26]. Here, we 
undertake a model-independent approach without speci- 
fying an underlying WIMP theory and using and a^j 
as our phenomenological parameters - for this reason we 
shall not address hadronic uncertainties (hidden in cr^j). 
A comprehensive work complementary to ours and done 
in the supersymmetric framework has been presented re- 
cently EaUH]. 

The paper is organised as follows. In the next section, 
we give some basic formulae for WIMP-nucleus recoil 
rates in direct detection experiments. In Section [III| the 
upcoming experimental capabilities are detailed, while 
Section [Tv| describes our Bayesian approach. We outline 
the relevant Galactic model uncertainties and our mod- 
elling of the velocity distribution function in Section [V] 
and present our results in Section VI before concluding 
in Section 



II. 



BASICS OF DIRECT DARK MATTER 
DETECTION 



Several thorough reviews on direct dark matter 
searches exist in the literature [1-6 . In this section, we 
simply recall the relevant formulae, emphasizing the im- 
pact of target properties and unknown quantities. 

The elastic recoil spectrum produced by WIMPs of 
mass and local density po on target nuclei N{A, Z) 
of mass rriN is 



dR 



(Er) 



Po 



(1) 

where v is the WIMP velocity in the detector rest frame, 
Ve is the Earth velocity in the Galactic rest frame, f{w) is 



the WIMP velocity distribution in the galactic rest frame 
and cTx-N is the WIMP-nucleus cross-section. The inte- 
gral is performed over V : v > Vmin(.Efj), where Vmin 
is the minimum WIMP velocity that produces a nuclear 
recoil of energy Er. Eq. ([T]) simply states that the recoil 
rate is the flux of WIMPs /^ow/"^x' averaged over the ve- 
locity distribution f{w), times the probability of interac- 
tion with one target nucleus a^-N- Anticipating the scale 
of future detectors, we will think of measuring dR/dEn in 
units of counts/ton/yr/keV. For non-relativistic (elastic) 
collisions - as appropriate for halo WIMPs, presenting 
v/c ^ 10^^ - the kinematics fixes the recoil energy 



ER{mx,v,A,e') = 
and the minimum velocity 



rriM 



{mx,ER,A) 



I uinEr 



in which 9' is the scattering angle in the centre of mass 



and fiN 



is the WIMP-nucleus reduced mass. 



In principle, all WIMP-nucleus couplings enter in the 
cross-section a^-N- However, we shall focus solely on 
spin-independent (SI) scalar interactions so that 



dEn 



niN 



2fi 



a'^jF\A,En) 



where a^j — 



^ [Zfp + {A- Z)Uf is the WIMP- 
nucleus spin-independent cross-section at zero momen- 
tum transfer and F{A, En) is the so-called form factor 
that accounts for the exchange of momentum. Assuming 
that the WIMP couplings to protons and neutrons are 

similar, /p ~ /„, and defining a^j = *^f^ being the 
WIMP-proton reduced mass), one gets 



da 



X-N 



dEr- 



<jPjA'F\A,En) 



(2) 



For the form factor, we use the parameterisation in [2] 
appropriate for spin-independent couplings, namely 

irt A TP \ qSin(g7-n) " (9fn)cos(gr„) ^ 
F{AEr) = 3 {qj^ exp(-(gs) /2) , 

with qr = 6.92 X IQ-^A^^^iE/keVy/^r/fm, s ~ 0.9 fm, 
rl = c2 + |7r2a2-5s2, c/fm = 1.23/11/3-0.6 and a ~ 0.52 
fm. 

As noticed above, in Eq. ([T]) is the Earth velocity 
with respect to the galactic rest frame and amounts to 
Ve = visr + Vpec + Vorb, whcrc Visr ~ 0(250) km/s is the 
local circular velocity, Vpec ~ 0(10) km/s is the peculiar 
velocity of the Sun (with respect to visr) and Vorb ~ 
0(30) km/s is the Earth velocity with respect to the Sun 
(i.e. the Earth orbit). Here, we are not interested in 
the annual modulation signal nor directional signatures 
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but rather in the average recoil rate - therefore we shall 
neglect Vpec and Vorb and take Ve — visr = const. 

Under these assumptions, Eq. ([ij may be recast in a 
very convenient way: 

{v7nin{m^,ER,A),Ve;Vo,Vesc) , (3) 

where we have used Eq. ([2|, defined 
Jv 

and made explicit the dependence of T on the velocity 
distribution parameters and Vesc- Below we discuss 
in more detail the connection between the parameters wq 
and visr- The distribution of DM is encoded in the fac- 
tor J- (and po)? whereas the detector-related quantities 
appear in A^F'^{A,Eii) (and Umm). The apparent de- 
generacy along the direction po<7gj/my. = const may be 
broken by using different recoil energies and/or different 
targets since J-" is sensitive to a non-trivial combination 
of iriy,, Er and A. Nevertheless, for very massive WIMPs 
^ rriN ^ O(IOO) GeV 3> rup, the minimum velocity 
becomes independent of m^, Vmin — \/ Eji/ (2m7v), and 
the degeneracy pQU^gj/m^ cannot be broken. Depending 
on the target being used, this usually happens for WIMP 
masses above a few hundred GeV. 

Ultimately, the observable we will be interested in is 
the number of recoil events in a given energy bin Ei < 
Er < E2: 

f^^ dR 
NRiE,,E2)^ dEne^ff— , (5) 

Eeff being the effective exposure (usually expressed in 
tonxyr) and dR/dER the recoil rate smeared according 
to the energy resolution of the detector ^{E), 

dR [ dR ( {E-E'f \ 

dE-R-j'^'dE-R^''^ v^^r"" 

Three fiducial WIMP models will be used to assess 
the capabilities of future direct detection experiments: 
=25, 50 and 250 GeV, all with a^j = 10"^ pb. These 
models are representative of well-motivated candidates 
such as neutralinos in supersymmetric theories |29j . 

III. UPCOMING EXPERIMENTAL 
CAPABILITIES 

Currently, the most stringent constraints on the SI 
WIMP-nucleon coupling are those obtained by the 
CDMS |30| and XENON [9 collaborations. While 
XENONIOO should probe the cross-section region down 
to 5 X 10"'*^ cm^ with data already in hand, the 



XENONIT [3T] detector, whose construction is sched- 
uled to start by mid 2011, is expected to reach another 
order of magnitude in sensitivity improvement. To test 
the agj region down to 10^"'^ cm^ = 10~^^ pb and below, 
a new generation of detectors with larger WIMP target 
masses and ultra-low backgrounds is needed. Since we 
are interested in the prospects for detection in the next 5 
to 10 years, we discuss new projects that can realistically 
be built on this time scale, adopting the most promising 
detection techniques, namely noble liquid time projec- 
tion chambers (TPCs) and cryogenic detectors operated 
at niK temperatures. 

In Europe, two large consortia, DARWIN and EU- 
RECA 133], gathering the expertise of several groups 
working on existing DM experiments are funded for R&D 
and design studies to push noble liquid and cryogenic ex- 
periments to the multi-ton and ton scale, respectively. 
DARWIN is devoted to noble liquids, having as main 
goal the construction of a multi-ton liquid Xe (LXe) 
and/or liquid Ar (LAr) instrument [M], with data tak- 
ing to start around 2016. The XENON, ArDM and 
WARP collaborations participate actively in the DAR- 
WIN project. EURECA is a design study dedicated to 
cryogenic dark matter detectors operated at mK tem- 
peratures. The proposed roadmap is to improve upon 
CRESST [3S] and EDELWEISS [35] technologies and 
build a ton-scale detector by 2018, with a SI sensitivity of 
about 10""*^ cm^ = 10^^° pb. The complementarity be- 
tween DARWIN and EURECA is of utmost importance 
for dark matter direct searches since a solid, uncontro- 
versial discovery requires signals in distinct targets and 
preferentially distinct technologies. In an international 
context, two engineering studies (MAX [37] and LZS [38J 
are funded in the US for ton to multi-ton scale LXe and 
LAr TPCs and the SuperCDMS/GEODM collaboration 
[35] plans to operate an 1.5 ton Ge cryogenic experiment 
at DUSEL W. In Japan, the XMASS experiment [H], 
using a total of 800 kg of liquid xenon in a single-phase 
detector, is under commissioning at the Kamioka under- 
ground laboratory [42 , while a large single-phase liquid 
argon detector, DEAP-3600 [43j, using 3.6 tons of LAr is 
under construction at SNOLab [44] . 

Given these developments, we will focus on the three 
most promising targets: Xe and Ar as examples of noble 
liquid detectors, and Ge as a case-study for the cryogenic 
technique. In the case of a Ge target, we assume an 
1.5 ton detector (Iton as fiducial target mass), 3 years 
of operation, an energy threshold for nuclear recoils of 
Ethr.Ge = 10 kcV and an energy resolution given by 

aceiE) = ViO.S)^ + (0.06)2£;/keV keV . (6) 

For a liquid Xe detector, we assume a total mass of 
8 tons (5 tons in the fiducial region), 1 year of operation, 
an energy threshold for nuclear recoils of Ethr.Xe — 10 
keV and an energy resolution of 

axe{E) = 0.6 keVV-B/keV . (7) 
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target 


e [tonxyrj 


r]cut 




eeff [tonxyr] Ethr |keVJ 


a{E) [ke 


VJ background events/ e^ff 


Xe 
Ge 
Ar 


5.0 
3.0 
10.0 


0.8 
0.8 
0.8 


0.5 
0.9 
0.8 


2.00 10 
2.16 10 
6.40 30 


Eq. ( 
Eq. ( 
Eq. ( 


7 
B 
S 


1 <1 
1 <1 
1 < 1 



TABLE I: Characteristics of future direct dark matter experiments using xenon, germanium and argon as target nuclei. In all 
cases the level of background in the fiducial mass region is negligible for the corresponding effective exposure. See Section |III| 
for further details. 



Finally, for a liquid Ar detector, we assume a total 
mass of 20 tons (10 tons in the fiducial region), 1 year 
of operation, an energy threshold for nuclear recoils of 
Ethr,Ar — 30 keV and an energy resolution of [45^ 



aAr{E) = 0.7 keVy/E/keV 



(8) 



To calculate realistic exposures, we make the following 
assumptions: nuclear recoils acceptances Ajyn of 90%, 
80% and 50% for Ge, Ar and Xe, respectively, and an 
additional, overall cut efficiency rjcut of 80% in all cases, 
which for simplicity we consider to be constant in energy. 
We hypothesise less than one background event per given 
effective exposure fe//, which amounts to 2.16 tonxyr in 
Ge, 6.4 tonxyr in Ar and 2 tonxyr in Xe, after allow- 
ing for all cuts. Such an ultra-low background will be 
achieved by a combination of background rejection using 
the ratio of charge-to-light in Ar and Xe, and charge-to- 
phonon in Ge, the timing characteristics of raw signals, 
the self-shielding properties and extreme radio-purity of 
detector materials, as well as minimisation of exposure 
to cosmic rays above ground. 

The described characteristics are summarised in Table 
|lj We note that in the following we shall consider recoil 
energies below 100 keV only; to increase this maximal 
value may add some information but the effect is likely 
small given the exponential nature of WIMP-induced re- 
coiling spectra. 



IV. STATISTICAL METHODOLOGY 

We take a Bayesian approach to parameter inference. 
We begin by briefly summarizing the basics, and we refer 
the reader to [IS] for further details. Bayesian inference 
rests on Bayes theorem, which reads 



pie\d) 



Pid\e)p{e) 
Pid) 



where p{Q\d) is the posterior probability density func- 
tion (pdf) for the parameters of interest, 0, given data 
d, p{d\Q) — £(6) is the likelihood function (when viewed 
as a function of 8 for fixed data d) and p{0) is the prior. 
Bayes theorem thus updates our prior knowledge about 
the parameters to the posterior by accounting for the in- 
formation contained in the likelihood. The normalization 
constant on the r.h.s. of Eq. ^ is the Bayesian evidence 
and it is given by the average likelihood under the prior: 



The evidence is the central quantity for Bayesian model 
comparison [4 7) , but it is just a normalisation constant 
in the context of the present paper. 

The parameter set 8 contains the DM quantities we 
are interested in (mass and scattering cross-section) , and 
also the Galactic model parameters, which we regard as 
nuisance parameters, entering the calculation of direct 
detection signals, namely po, vq, Vesc, k, see Eq. ([3| and 
Sectionjv] We further need to define priors p(8) for all of 
our parameters. For the DM parameters, we adopt flat 
priors on the log of the mass and cross-section, reflecting 
ignorance on their scale. For the Galactic model param- 
eters, we choose priors that reflect our state of knowl- 
edge about their plausible values, as specified in the next 
section. Those priors are informed by available observa- 
tional constraints as well as plausible estimations of un- 
derlying systematical errors, for example for po- Finally, 
the likelihood function for each of the direct detection ex- 
periments is given by a product of independent Poisson 
likelihoods over the energy bins: 



(11) 



where Nb is the number of counts in each bin (generated 
from the true model with no shot noise, as explained be- 
low) and Nr = A^fl(£;™'", S™'''') is the number of counts 
in the 6-th bin (in the energy range i?™'" < E < E^^^) 
when the parameters take on the value 8, and it is given 
byEq. We use 10 bins for each experiment, uniformly 
spaced on a linear scale between the threshold energy and 
100 keV. We have checked that our results are robust if 
we double the number of assumed energy bins. Using the 
experimental capabilities outlined in Section III we com- 



p{d) = / dQp{d\Q)p{@). 



(10) 



pute the counts Nr that the benchmark WIMPs would 
generate, and include no background events since the ex- 
pected background level in the fiducial mass region is 
negligible (cf . Table |l]) . The mock counts are generated 
from the true model, i.e. without Poisson scatter. This 
is because we want to test the reconstruction capabilities 
without having to worry about realization noise (such a 
data set has been called "Asimov data" in the particle 
physics context [48j ) . 

To sample the posterior distribution we employ the 
MultiNest code [15H51| . an extremely efficient sampler 
of the posterior distribution even for likelihood functions 
defined over a parameter space of large dimensionality 
with a very complex structure. In our case, the likeli- 
hood function is unimodal and well-behaved, so Monte 
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Parameter 


Prior range Prior constraint 


logio ((^Sj/Pb) 


(0.1,3.0) Uniform prior 
(—10,-6) Uniform prior 


Po/(GeV/cm") 
Wo/ (km/s) 
«esc/ (km/s) 
k 


(0.001,0.9) Gaussian: 0.4 ±0.1 
(80,380) Gaussian: 230 ± 30 
(379,709) Gaussian: 544 ± 33 
(0.5,3.5) Uniform prior 



TABLE II: Paramete rs us ed in our analysis, with their prior range (middle column) and the prior constraint adopted (rightmost 
column). See Section IV and Iv] for further details. 



Carlo Markov Chain (MCMC) techniques would be suf- 
ficient to explore it. However, MultiNest also computes 
the Bayesian evidence (which MCMC methods do not re- 
turn), as it is an implementation of the nested sampling 
algorithm (55] . In this work, we run MultiNest with 2000 
live points, an efficiency parameter of 1.0 and a tolerance 
of 0.8 (see [491E0] for details). 



V. VELOCITY DISTRIBUTION AND 
GALACTIC MODEL PARAMETERS 

We now move onto discussing our modeling of the ve- 
locity distribution function and the Galactic model pa- 
rameters that are input for Eq. ([s]). We model only 
the smooth component of the velocity distribution - 
recent results from numerical simulations indicate that 
the velocity distribution component arising from lo- 
calised streams and substructures is likely sub-dominant 
in the calculation of direct dark matter detection sig- 
nals [531 El]. 

We model the velocity distribution function as spheri- 
cal and isotropic, and parameterise it as |55j . 



exp 



\ kvl ) 



- 1 



if W < Vesc ("22) 
if W > Vesc 



This velocity distribution function was found to be flex- 
ible enough to describe the range of dark matter halo 
profiles found in cosmological simulations [SS] . Boosting 
into the rest frame of the Earth implies the transforma- 
tion — u-^+Wg +2ti?;eCos0, where is the angle between 
V and Ve ~ visr- The shape parameter that determines 
the power law tail of the velocity distribution is k, the 
escape velocity is t^esc, while vq is a fit parameter that we 
discuss in detail below, and Nj is the appropriate nor- 
malisation constant. The special case k — 1 represents 
the standard halo model with a truncated Maxwellian 
distribution, and the corresponding expressions for Nf 
and J- have been derived analytically in the literature - 
see for instance [T51. Note as well that, for any value 
of fc, this distribution matches a Maxwellian distribution 
for sufficiently small velocities w and if Vesc > vo- 

The high- velocity tail of the distributions found in nu- 
merical simulations of pure dark matter galactic halos are 
well modelled by 1.5 < k < 3.5 ^5]- our analysis we 
will expand this range to also include models that behave 



similar to pure Maxwellian distributions near the tail of 
the distribution, so that in our analysis we vary k in the 
range 



fc = 0.5 - 3.5 (flat) 



(13) 



We adopt an uniform (i.e., flat) prior within the above 
range for k. 

The range we take for the Vesc is motivated by the re- 
sults of Ref . [56J , where a sample of high- velocity stars is 
used to derive a median likelihood local escape velocity 
of Vesc — 544 km/s and a 90% confidence level interval 
498 km/s < v^sc < 608 km/s. Assuming Gaussian errors 
this translates into a la uncertainty of 33 km/s. It is im- 
portant to note that this constraint on the escape velocity 
is derived assuming a range in the power law tail for the 
distribution of stars in the local neighbourhood, which 
is then related to the power law tail in the dark matter 
distribution [56] . Motivated by obtaining conservative 
limits on the reconstructed mass and cross-section of the 
dark matter, in our modelling we will not include such 
correlations between the escape velocity and the power 
law index fc, so that in the end we take a Gaussian prior 
on Vesc with mean and standard deviation given by 



544 ±33 km/s (Icr) 



(14) 



Having specified ranges for Vesc and k, it remains to 
consider a range for vq in Eq. (12). As defined in that 
equation, the quantity vq does not directly correspond 
to the local circular velocity, visr, but rather is primarily 
set by Visr and the dark matter profile. Following a pro- 
cedure similar to that discussed in Ref. [5S], we find the 
range of values vq compatible with a given a dark matter 
halo profile, po and a range for visr- For the above range 
in Visr and the values po in Eq. (16) below, we find that 
the parameter vq can take values in the range 200 — 300 
km/s for pure Navarro-Frenk- White (NEW) dark matter 
halos with outer density slopes p cx . Larger values of 
Wo are allowed for steeper outer density slopes, though the 
range is found to not expand significantly if we restrict 
ourselves to models with outer slopes similar to the NEW 
case. With these caveats in mind regarding the mapping 
between vq and visr for steeper outer slopes, for simplic- 
ity and transparency in our analysis, we will consider a 
similar range for vq as for the local circular velocity, so 
we take vq = visr (that holds in the case of the standard 
halo model). 
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FIG. 1: The joint 68% and 95% posterior probability contours in the — a^j plane for the three DM benchmarks 
(m^ = 25, 50, 250 GeV) with fixed Galactic model, i.e. fixed astrophysical parameters. In the left frame we show the re- 
construction capabilities of Xe, Ge and Ar configurations separately, whereas in the right frame the combined data sets Xe+Ge 
and Xe+Ge+Ar are shown. 



For the local circular velocity and its uncertainty, a va- 
riety of measurements presents a broad range of central 
values and uncertainties [57] . To again remain conserva- 
tive we use an interval bracketing recent determinations: 



Vo = Vis 



230 ±30 km/s (la) 



(15) 



where we take a Gaussian prior with the above mean and 
standard deviation. To account for the variation of the 
local density of dark matter in our modeling, we will take 
a mean value and error given by |58|, I59j 



(16) 



po = 0.4±0.1 GeV/cm'' (Icr) 



There are several other recent results that determine po, 
both consistent [6Uj and somewhat discrepant |61| with 
our adopted value. Even in light of these uncertainties, 
we take Eq. ( 16 1 to represent a conservative range for the 



purposes of our study. 

For completeness Table [TT| summarises the information 
on the parameters used in our analysis. 



VI. RESULTS 

A. Complementarity of targets 

We start by assuming the three dark matter bench- 
mark models described in Section |ll] (m^ — 25, 50, 250 
GeV with a^j — 10^^ pb) and fix the Galactic model 
parameters to their fiducial values, po = 0.4 GeV/cm'^, 
Vq = 230 km/s, v^g^ — 544 km/s, k = 1. With the exper- 
imental capabilities outlined in Section |III[ we generate 



mock data that in turn are used to reconstruct the poste- 
rior for the DM parameters and cr^j. The left frame 
of Fig. [T] presents the results for the three benchmarks 
and for Xe, Ge and Ar separately. Contours in the figure 
delimit regions of joint 68% and 95% posterior probabil- 
ity. Several comments are in order here. First, it is ev- 
ident that the Ar configuration is less constraining than 
Xc or Ge ones, which can be traced back to its smaller A 
and larger Ethr- Moreover, it is also apparent that, while 
Ge is the most effective target for the benchmarks with 
= 25, 250 GeV, Xe appears the best for a WIMP with 
= 50 GeV (see below for a detailed discussion). Let 
us stress as well that the 250 GeV WIMP proves very 
difficult to constrain in terms of mass and cross-section 
due to the high-mass degeneracy explained in Section [ll] 
Taking into account the differences in adopted values and 
procedures, our results are in qualitative agreement with 
Ref. |27| , where a study on the supersymmetrical frame- 
work was performed. However, it is worth noticing that 
the contours in Ref. [27] do not extend to high masses 
as ours for the 250 GeV benchmark - this is likely be- 
cause the volume at high masses in a supersymmetrical 
parameter space is small. 

In the right frame of Fig. [T]we show the reconstruction 
capabilities attained if one combines Xe and Ge data, 
or Xe, Ge and Ar together, again for when the Galac- 
tic model parameters are kept fixed. In this case, for 
= 25, 50 GeV, the configuration Xe-|-Ar-|-Ge allows 
the extraction of the correct mass to better than 0(10) 
GeV accuracy. For reference, the (marginalised) mass 
accuracy for different mock data sets is listed in Table 



III For = 250 GeV, it is only possible to obtain a 



lower limit on 
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FIG. 2: The joint 68% and 95% posterior probability contours in the — a^j plane for the case in which astrophysical 
uncertainties are taken into account. In the left frame, the effect of marginaUsing over po, vo and all four (po, vo, Vesc, k) 
astrophysical parameters is displayed for a Xe detector and the 50 GeV benchmark WIMP. In the right frame, the combined 
data sets Xe+Ge and Xe+Ge+Ar are used for the three DM benchmarks {rn^ — 25, 50, 250 GeV). 





Percent Ic 
= 25 GeV 


r accuracy 

= 50 GeV 


Xe 
Ge 
Ar 
Xe+Ge 
Xe+Ge+Ar 


6.5% (14.3%) 
5.5% (16.0%) 
12.3% (23.4%) 
3.9% (10.9%) 
3.6% (9.0%) 


8.1% (20.4%) 
7.0% (29.6%) 
14.7% (86.5%) 
5.2% (15.2%) 
4.5% (10.7%) 



TABLE III: Marginalised percent la accuracy of the DM mass reconstruction for the benchmarks — 25, 50 GeV. Figures 
between brackets refer to scans where the astrophysical parameters were marginalised over (with priors as in Table [ll| , while 
the other figures refer to scans with the fiducial astrophysical setup. 



Fig. [2] shows the results of a more reaUstic analysis, 
that keeps into account the large uncertainties associated 
with Galactic model parameters, as discussed in Section 
|Vj The left frame of Fig. [2] shows the effect of varying 
only pq (dashed lines, blue surfaces), only vq (solid lines, 
red surfaces) and all Galactic model parameters (dotted 
Hnes, yellow surfaces) for Xe and = 50 GeV. The 
Galactic model uncertainties are dominated by po and 
fo, and, once marginalised over, they blow up the con- 
straints obtained with fixed Galactic model parameters. 
This amounts to a very significant degradation of mass 
(cf. Tabic III I and scattering cross-section reconstruction. 
Inevitably, the complementarity between different targets 
is affected - see the right frame of Fig. [2j Still, for the 
50 GeV benchmark, combining Xe, Ge and Ar data im- 
proves the mass reconstruction accuracy with respect to 
the Xe only case, essentially by constraining the high- 
mass tail. 

In order to be more quantitative in assessing the use- 
fulness of different targets and their complementarity, we 
use as figure of merit the inverse area enclosed by the 
95% marginalised contour in the logiQ^m^) — logiQ^a^j) 
plane inside the prior range. Notice that for the 250 



GeV benchmark the degeneracy between mass and cross- 
section is not broken - this does not lead to a van- 
ishing figure of merit (i.e. infinite area under the con- 
tour) because we are restricting ourselves to the prior 
range. Fig. [3] displays this figure of merit for several 
cases, where we have normalised to the Ar target at 
= 250 GeV with fixed Galactic model. Analyses 
with fixed Galactic model parameters are represented by 
empty bars, while the cases where all Galactic model pa- 
rameters are marginalised over with priors as in Table [Tl| 
are represented by filled bars. Firstly, one can see that all 
three targets perform better for WIMP masses around 50 
GeV than 25 or 250 GeV if the Galactic model is fixed. 
When astrophysical uncertainties are marginalised over, 
the constraining power of the experiments becomes very 
similar for benchmark WIMP masses of 25 and 50 GeV. 
Secondly, Fig. J3] also confirms what was already appar- 
ent from Fig. [l Ge is the best target for = 25, 250 
GeV (althougnby a narrow margin) , whereas Xe appears 
the most effective for a 50 GeV WIMP (again, by a nar- 
row margin). Furthermore, the inclusion of uncertainties 
drastically reduces the amount of information one can 
extract from the data: the filled bars are systematically 
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FIG. 3: Figure of merit quantifying the relative information 
gain on Dark Matter parameters for different targets and com- 
binations thereof. The values of the figure of merit are nor- 
malised to the Ar case at — 250 GeV with fixed astrophys- 
ical parameters. Empty (filled) bars are for fixed astrophysical 
parameters (including astrophysical uncertainties). 



below the empty ones. Now, astrophysical uncertainties 
affect the complementarity between different targets in a 
non-trivial way. To understand this point, let us focus 
on the two rightmost bars for each benchmark in Fig. |3j 
corresponding to the data sets Xe-j-Ge and Xe-|-Ge-|-Ar. 
For instance, in the case of a 250 GeV WIMP, astrophysi- 
cal uncertainties seem to reduce target complementarity: 
adding Ar to Xe-fGe leads to a significant increase in 
the figure of merit for analyses with fixed astrophysics 
(empty bars) but has a negligible effect for analyses with 
varying astrophysical parameters (filled bars). For low 
mass benchmarks, the effect of combining two (Xe-fGe) 
or three targets (Xe-|-Ge-|-Ar) is to increase the figure of 
merit by about a factor of 2 compared to Xe alone or Ge 
alone, almost independently of whether the astrophysical 
parameters are fixed or marginalised over. However, the 
overall information gain on the Dark Matter parameters 
(for light WIMPs) is reduced by a factor ~ 10 if astro- 
physical uncertainties are taken into account, compared 
to the case where the Galactic model is fixed. 



B. Reduction in uncertainties and self-calibration 

The uncertainties used thus far and outlined in Section 
[V| are a reasonable representation of the current knowl- 
edge. For illustration it is also interesting to consider 
the effect of tighter constraints on Galactic model pa- 
rameters in the reconstruction of WIMP properties. We 
start by computing the correlation coefficient between 
the parameters (m^, a^j, po, Vq, v^sc^ k) when they 
are constrained by the combined data set Xe-|-Ge-|-Ar 



FIG. 4: The effect of reducing the uncertainty on the astro- 
physical parameters po and vq. The red surfaces refer to the 
scan using the fiducial astrophysical setup; the yellow sur- 
faces (and dotted lines) indicate the effect of marginalising 
over the uncertainties in Table [ll] the blue surfaces (and solid 
lines) correspond to the reduced uncertainties po = 0.4 ± 
0.028 GeV/cm^ vq = 230±9.76 km/s, v^sc = 544±33 km/s, 
fc = 0.5-3.5. 



- see Table |IV[ Glearly, for all benchmark models, a^j 
and Po as well as and vq are strongly anti-correlated. 
The anti-correlation between a^j and po is obvious since 
dR/dEj^ cx (TgjPQ- As for the degeneracy between 
and Vq, it is easy to verify that, for Vmin <C wo ^ 
Vpfir, T defined in Eq. (|4]) goes app roxi mately as I/wq and 
' ~ " also shows a small 



IV 



thus dR/dEji (X l/(m^wo). Table 
(anti-)correlation between cr^j and wo; all other correla- 
tions are negligible. Therefore, po and are the domi- 
nant sources of uncertainty and their more accurate de- 
termination will lead to a significant improvement on the 
reconstruction of m^^ and cr^j. To illustrate this point we 
follow [5^ and apply a 7% (4.2%) uncertainty on po (^^o), 
while maintaining the same central values as before, thus 
reducing the realistic error bars used above by a factor 
^ 3.0 — 3.5 for both parameters. The results are shown in 
Fig. [4] where we consider the combination Xe-j-Ge+Ar. A 
future more constrained astrophysical setup may indeed 
lead to a better reconstruction of the WIMP mass and 
scattering cross-section. 

To this point we have studied the impact of Galac- 
tic model uncertainties on the extraction of DM proper- 
ties from direct detection data. However, once a positive 
signal is well-established, it may be used to determine 
some of the Galactic parameters directly from direct de- 
tection data (see e.g. [H]), without relying on external 
priors. This would amount to achieving a self-calibration 
of the astrophysical uncertainties affecting direct detec- 
tion rates. In order to explore such possibility we re-ran 
our analysis but dropping the Gaussian priors on poj 
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-0.006 
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-0.004 




0.874 


-0.011 
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0.022 


V 

-SI 
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-0.887 


-0.237 0.116 


0.010 






-0.957 


-0.175 0.026 


-0.031 






-0.452 


-0.525 -0.024 


0.015 






0.013 -0.005 


0.005 








0.014 -0.010 


0.030 








0.002 0.015 


0.010 


VO 






- -0.087 


-0.004 








- -0.151 


0.011 








- -0.049 


-0.008 










0.000 










-0.009 










0.001 



TABLE IV: The correlation factors r{X,Y) — cov{X, Y)/{a{X)a{Y)) for the posteriors obtained from the combined data set 
Xe-|-Ge-|-Ar and including the astrophysical uncertainties with priors as in Table [ill 
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FIG. 5: The marginalised posterior distribution function for (left frame) and vo (right frame) with the data sets Xe, Xe-j-Ge 
and Xe+Ge-f-Ar for the 50 GeV benchmark. The parameters po, vo, v^sc and k were varied in the ranges indicated in the 
middle column of Table [H] with a uniform prior and no constraint on astrophysics was applied. The probability distributions 
are therefore a result of trie constraining power of direct detection data only, which have the potential to achieve self-calibration 
of the circular velocity. 



and Wesc described in Section |V] Instead, we used uni- 
form, non- informative priors on pQ, vq, Vesc and k in the 
ranges indicated in the middle column of Table We 
focus on the 50 GeV benchmark and use the data sets 
Xe, Xe-|-Ge and Xe-f Ge-j-Ar. With this large freedom on 
the astrophysical side, it turns out that direct detection 
data alone leave poi Vesc and k unconstrained within their 
ranges while a^j is pinpointed within approximately one 
order of magnitude. Only the DM mass and the cir- 
cular velocity vq can be constrained by direct detection, 
as shown in Fig. [5] This figure stresses two interest- 
ing results. First, if ~ 50 GeV (and a^j — 10~^ 
pb), the next generation of experiments will be able to 
determine the WIMP mass within a few tens of GeV 
(percent la accuracy of 11.8%) even with very loose as- 
sumptions on the local DM distribution. Second, the 
right frame in Fig. [5] shows that the combination of Xe, 
Ge and Ar targets is very powerful in constraining vq on 
its own without external priors. In particular, the data 
set Xe-|-Ge-|-Ar (solid blue line) is sufficient to infer at 
1(7 Vq = 238 ± 22 km/s (compared to the top-hat prior 
in the range 80—380 km/s). This represents already a 
smaller uncertainty than the present-day constraint that 



we have taken, vq — 230 ± 30 km/s - in case of a positive 
signal, a combination of direct detection experiments will 
probe in an effective way the local circular velocity. Re- 
peating the same exercise for the 25 GeV benchmark we 
find good mass reconstruction but a weaker constraint: 
Vq — 253 ± 39 km/s. Again, we stress that the quoted vq 
uncertainties in this paragraph do not take into account 
possible systematic deviations from the parameterisation 
inEq. (pi). 



VII. CONCLUSIONS 

We have discussed the reconstruction of the key phe- 
nomenological parameters of WIMPs, namely mass and 
scattering cross-section off nuclei, in case of positive de- 
tection with one or more direct DM experiments planned 
for the next decade. We have in particular studied the 
complementarity of ton scale experiments with Xe, Ar 
and Ge targets, adopting experimental configurations 
that may realistically become available over this time 
scale. 

To quantify the degree of complementarity of different 
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targets we have introduced a figure of merit measuring 
the inverse of the area enclosed by the 95% marginahsed 
contours in the plane \og^Q{m^) — \ogiQ{agj). There is a 
high degree of complementarity of different targets: for 
our benchmark with = 50 GeV and our fiducial set of 
Galactic model parameters, the relative error on the re- 
constructed mass goes from 8.1% for an analysis based on 
a xenon experiment only, to 5.2% for a combined analysis 
with germanium, to 4.5% adding also argon. Allowing 
the parameters to vary within the observational uncer- 
tainties significantly degrades the reconstruction of the 
mass, increasing the relative error by up to a factor of ~4 
for xenon and germanium, especially due to the uncer- 
tainty on pq and vq. However, we found that combining 
data from Ar, Ge and Xe should allow to reconstruct a 50 
GeV WIMP mass to 11.8% accuracy even under weaker 
astrophysical constraints than currently available. 

Although the mass reconstruction accuracy may ap- 
pear modest, any improvement of this reconstruction is 
important, in particular in view of the possible measure- 
ment of the same quantity at the Large Hadron Collider 
at CERN. The existence of a particle with a mass com- 
patible, within the respective uncertainties, with that de- 
duced from direct detection experiments would provide a 
convincing proof that the particles produced in acceler- 
ators are stable over cosmological time scales. Although 
this is not sufficient to claim discovery of DM [25j, it 



would certainly be reassuring. 

Despite the strong dependence of direct detection ex- 
periments on the Galactic model degrades the reconstruc- 
tion of DM properties, it does open up the possibility to 
potentially constrain the local distribution of DM, in case 
of detection with multiple targets. For example in the 
case of a low mass 50 GeV WIMP, we have shown that 
the local circular velocity can be determined from direct 
detection data alone more accurately than it is presently 
measured using the local distribution of stars and gas 
clouds. Additionally, directly detecting DM provides the 
most realistic way of measuring the local DM velocity 
distribution. This will in principle provide invaluable in- 
formation on the structure and formation of the Milky 
Way halo. 
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